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SUMM  AR  Y 


In  preparation  for  a computer  simiilation  demonstration  of  non- 
stationary post-detection  compensation,  we  have  concerned  ourselves  in 
this  report  with  the  task  of  generating  a one -dimensional  image  signal 
with  a H-i  power  spectrum,  matching  a two-dimensional  image  with  a 
power  spectrum.  Using  the  insight  provided  by  our  previous  work 
on  nonstationary  post- detection  compensation,  we  formulate  the  image 
as  a set  of  rectangular  pulses  of  various  widths,  randomly  distributed 
across  the  field-of-view.  By  relating  the  distribution  of  pulse  height  to 
pulse  width,  we  are  able  to  keep  the  average  pulse  height  independent  of 
width,  but  make  the  second  moment  of  pulse  height  depend  on  pulse 
width.  This  allows  us  to  adjust  the  average  pulse  count  as  a function 
of  width  so  that  the  power  spectrum  is  , and  yet  maintain  that 

approximately  one-half  of  the  "object  region"  is  empty  — a general 
characteristic  of  satellite  images,  and  a key  item  in  exploitation  of 
nonstationary  post- detection  compensation. 
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Introduction 


In  a previous  report,  we  developed  the  concept  of  nonstationary 
post- detection  compensation  as  a way  to  take  account  of  the  information 
available  a posteriori  in  defining  a filter  for  compensating  an  image.  The 
nonstationary  filter  concept  as  presented  there  was  developed  by  means  of 
a combination  of  analytic  and  heuristic  techniques.  This  left  the  demon- 
stration of  the  concept's  suitability  as  a matter  for  further  study.  As  a 
practical  matter,  we  believe  that  the  only  really  sound  approach  to  this 
demonstration  is  an  empirical  one.  Accordingly,  we  have  sought  to  apply 
our  nonstationary  post-detection  compensation  concept  to  a one- dimensional 
problem,  so  as  to  facilitate  computer  demonstration. 


It  is  the  objective  of  this  report  to  define,  in  sufficient  detail  to 
support  subsequent  computer  programming,  a procedure  for  generating 
appropriate  one-dimensional  images  manifesting  appropriate  target  image 
statistics.  As  noted  earlier,  the  key  feature  of  a realistic  two-dimensional 
target  image  in  a region  where  an  interesting  abundance  of  detail  is  present 
is  the  fact  that  its  power -spectrum  behaves  as  h"®  , where  h denotes  a 

spatial  frequency,  and  the  power -spectrum  is  the  fully  two-dimensional 
power -spectrum.  In  restricting  ourselves  to  treatment  of  the  one -dimensional 
image  problem,  we  cause  the  power- spectrum  of  interest  to  take  on  a 

-dependence.  In  the  next  section,  we  shall  discuss  the  other  character- 
istics we  seek  in  the  image  pattern  we  shall  generate.  In  the  sections  after 
that,  we  shall  present  a method  of  generating  such  an  image,  and  then  will 
prove  that  the  one -dimensional  image  so  generated  does  indeed  have  the 
desired  properties. 


D.  L.  Fried,  "Analysis  of  Techniques  for  Imaging  Through  the  Atmos 
phere,"  Rome  Air  Development  Center  Report  No.  RADC-TR-77-196, 
June  1977;  Chapter  2,  AD#  B019  482L. 


2.  Desired  Image  Properties 

First  and  foremost,  the  one -dimensional  image  we  seek  to  gener- 
ate will  have  a power -spectrum  behaving  as  k"*-  . This  leads  us  to  hypo- 

thesize that  the  image  should  be  considered  to  be  formed  of  a randomly 
dispersed  set  of  rectangular  pulses,  with  the  number  of  pulses  of  each 
width  and  the  average  intensity  for  each  pulse  so  chosen  that  the  composite 
pattern  has  a power- spectrum.  We  shall  incorporate  this  hypothesized 

formulation  of  the  total  image  pattern  as  a sum  of  rectangular  pulses  into 
our  presentation  of  the  other  desired  image  properties. 


We  consider  the  object  pattern  to  be  formed  of  a randomly  selected 
array  of  pulses.  We  expect  the  pulse  sizes  to  cover  a range  from  some 
maximum  size,  comparable  to  about  one -eighth  of  the  total  image  field- 
of-view,  to  a minimum  size.  The  minimiun  size  is  chosen  simply  for 
lack  of  interest  in  smaller  details,  with  the  understanding,  however,  that 
there  actually  is  no  smallest  size  cut-off. 


We  assume  that  the  number  of  pulses  of  each  width  is  a random 
variable,  as  is  the  positioning  of  each  pulse  and  the  intensity  of  each 
pulse.  These  values  are,  however,  subject  to  certain  constraints.  To 
maintain  the  compactness  of  the  total  image  within  the  field-of-view,  i.  e.  , 
the  desired  property  that  the  image  appear  to  be  that  of  a contiguous  whole 
rather  than  of  a uniform  random  distribution  over  the  entire  available  field 
of-view,  we  assume  some  constraint  on  the  distribution  of  positioning  of 
the  pulses.  To  maintain  the  concept  of  the  target  object  as  made  up  of 
different  size  elements  of  related  characteristics,  we  shall  assume  that 
the  average  intensity  (i.e,  , the  pulse  height)  associated  with  each  pulse 
width  be  independent  of  the  pulse  width.  (However,  we  allow  the  distribu- 
tion of  pulse  height  to  be  a function  of  pulse  widths  — providing  that  the 
function  does  not  affect  the  average  pulse  height.  ) 
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The  average  number  of  pulses  of  each  width  is  constrained  by  the 
need  to  achieve  a power-spectrum,  but  to  the  extent  that  the  relation- 

ship between  the  mean  and  second  moment  of  the  pulse  intensity  can  be 
made  a function  of  the  pulse  width,  we  can  correspondingly  adjust  the  mean 
number  of  pulses  as  a function  of  pulse  width.  Our  concern,  in  remarking 
on  this  fact,  is  that  examination  of  typical  real  objects  of  interest  shows 
that  in  the  compact  region  which  contains  most  of  a two-dimensional  object, 
about  three-quarters  of  the  space  is  entirely  empty.  For  a one -dimensional 
object  pattern,  the  corresponding  factor  would  be  that  about  one-half  of  the 
space  will  be  entirely  empty.  In  our  choice  of  the  intensity  (i.e.  , pulse 
height)  distribution  for  each  element  size  (i.e.  , for  each  pulse  width),  we 
shall  choose  the  distribution  so  that  the  second  moment  implies  an  average 
number  of  elements  (i.e,  , pulses)  of  that  size,  which  when  "suftimed"  over 
all  element  sizes  considered  results  in  the  desired  50%  vacancy. 

It  should  be  remarked  that  this  procedure  only  specifies  the  average 
number  of  elements  of  each  size.  The  actual  number  of  elements  will  be  a 
random  variable  with  this  average  value  — except  that  for  the  largest  ele- 
ment, we  shall  impose  the  requirement  that  there  is  always  one,  and  only 
one  element.  This  will  not  be  a random  variable. 

With  the  image  properties  as  set  forth  here,  we  are  now  ready  to 
consider  definition  of  a procedure  for  generating  a randomly  "chosen" 
sample  image  with  these  properties.  This  is  taken  up  in  the  next  section. 


ri 
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3.  Sample  Image  Generation 

We  propose  to  consider  an  image  defined  on  a "field-of -view"  of 
2»  points  (where  in  practice  we  shall  use  n = 10  , 2"  = 1024).  We  num- 

•| 

, her  these  positions  as  0 . 1 , 2 , . . . 2»  - 1 . To  maintain  the  compact- 

ness of  the  object,  we  define  a standard  deviation,  , defining  the  image 
pattern  spread,  and  consider  the  center  point  of  each  pulse  making  up  the 
image  to  be  randomly  distributed  in  accordance  with  a gaussian  probability 
distribution  with  this  standard  deviation  and  mean  value  I =2<*-i)  . Thus, 
if  some  pulse  is  centered  at  i , then  I will  be  chosen  as  a gaussian  ran- 
dom variable  with  probability  density 

Pc  (i)  = {2tt  a,»F*  exp  [ (jl-7)a/a,3  ] (D 

We  will  use  or  for  n = 10  , Oj  = 2®  = 256,  and  T = 512  . This 

means  that  the  nominal  width  of  our  object  pattern  will  be  spread  over  about 
3 a,  = 768  points,  or  three-quarters  of  the  field-of -view,  and  that  the 
pattern  will  be  nominally  centered  in  the  middle  of  our  field-of-view. 

We  shall  consider  our  largest  element  to  be  a. pulse  of  width 

*o=2<*-«>=  27  = 128  (2) 

We  shall  consider  a set  of  pulse  widths 

a^  = , for  W=0,l,2,...,  n-3  , (3) 

so  that  for  the  n = 10  size  field-of-view  we  intend  to  consider,  w takes 
8 values  (i.  e.  , W = 0 , 1,2 7 ) . 

To  interpret  the  requirement  that  approximately  one-half  of  the 
object  area  be  entirely  empty  of  pulse  elements,  we  need  first  of  all  to 
define  what  we  mean  by  the  object  area.  Arguing  from  insight  gained  from 
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a study  of  sample  targets  with  K~<t  power  spectra,  we  would  suggest  that 
the  object  area  be  considered  to  be  about  six  times  the  size  of  the  largest 
element.  Within  this  size  region,  the  target  elements  take  up  about  one- 
half  of  the  space  in  one -dimension  (and  one -quarter  of  the  space  in  two- 
dimensions).  This  means  that  the  object  area  is  6 a^  . Assuming  a power 
law  dependence  for  the  average  number  of  elements  of  size  , the 

average  number  is 

= bH  , (4) 

then  the  average  space  taken  up  by  all  the  elements  will  be 

•''i**"-  'I 

H«0  MsO 

7 

= ^ sy,  (b/2)M  (5)* 

NsO 

But  for  half  the  object  space,  i.  e. , half  of  6 a^  , to  correspond  to  ^ , 
it  follows  that 


3 = ^ (h/2f 

w«o 


(b/Z)"-! 
(b/2)  -1 


i 

l-(b/2) 


(6) 


* This  formulation  is  strictly  valid  only  where  the  area  density  of  elements 
is  so  low  that  tiiere  is  a negligible  probability  of  element  overlap.  With 
our  5C%  area  density  objective,  this  approximation  is  not  strictly  valid. 
However,  at  this  density  the  overlap  correction  factor  will  be  close  enough 
to  unity  that  we  can  ignore  it  — doing  so  being  no  greater  an  approxima- 
tion than  our  "estimation"  that  the  area  fill  factor  should  be  50%.  The  re- 
sult of  this  approximation  in  Eq.  (5)  will  be  that  the  patterns  we  generate 

will  actually  have  a slightly  less  than  50%  area  fill  factor. 


Solving  Eq.  (6)  for  b , we  get 


b = 1.  3333  . (7)  j 

\ 

which  implies  that  = { 1 , 1.3333  , 1.  7778  , 2.3704  , 3.1605  , 4.2140  , | 

5.6187  . 7.4915)  . j 

We  shall  show  in  our  analysis  that  for  this  set  of  values  of  to  i 

be  compatible  with  a h-^  one-dimensional  power-spectrum,  it  is  necessary  | 

p 

that  the  ratio  of  the  second  moment  of  pulse  height  (i.  e.  , of  element  inten-  v 

sity)  divided  by  the  square  of  the  first  moment,  i.  e.  , )*  » depend  , ■ 

on  W according  to  the  relationship 

= c-  . (8) 

where 

c = 1.  5 (9) 

A proof  of  this  will  be  presented  in  the  next  section.  In  order  to  achieve 
this  ratio  of  second  moment  to  first  moment  squared,  we  shall  assume 
that  the  pulse  heights  are  random  variables,  distributed  for  each  pulse 
width  according  to  a log-normal  distribution,  with  the  basic  parameters 
of  the  distribution  dependent  on  the  pulse  width. 

i 

The  log-normal  distribution  is  characterized  by  two  parameters.  ' 

The  first  is  the  logarithmic  standard  deviation,  , and  the  second  is  ' 

the  logarithmic  mean,  L . If  obeys  a log-normal  distribution  with 
these  parameters,  then  the  probability  density  associated  with  1^  is 

P,  (I,)  = {2tt  1^-*  exp  - L]a/o,*)  . (10) 

where  ^ is  a reference  intensity  value  corresponding  to  the  median 
value  of  intensity. 


<(iL> 

<In>» 
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Accordingly,  we  can  write  for  the  first  moment 


<!,>  = J dip, (I)  I 


= Ip  J di  (2n  exp 


U)  . 


where  we  have  introduced  tiie  notation  l and  replaced  the  variable 

of  integration  I by  £ . Completing  the  square  in  the  exponential,  this  can 


be  rewritten  as 


<I„>  = (2tt  o^a)-i/a  {•  dl  exp  r-^ ; 


(a^a  + Lya  . La 

XexpffJ  — ; J 


t '-P  |_‘i  ] 


It  follows  from  Eq.  (12)  that  the  first  moment,  (I„)  will  Ise  independent 
of  W (and  equal  to  the  mean  pulse  height  for  the  largest  element),  if 

L = a, a 

In  the  same  manner,  we  can  write  for  the  second  moment 

<!,» ) = f di  p,  a)  I* 

= l^a  f df  (2tt  a,aH/»  e*q>  ] exp  (2i)  . 


- 7 - 
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Completing  the  aquare  in  the  exponential,  this  can  be  rewritten  as 


, ^ U-2a»-E)*^ 

<!•)  = J dl  .xp  [-i 


« •‘P  [i  —^5 ] 


Combining  Eq.  's  (12)  and  (15),  and  making  use  of  Eq.  (13),  we  get 

^■L,)  - exp  la  * ) 

^s  «*P  '''i  • 

It  fellows  from  this  that  Eq.  's  (8)  and  (9)  will  be  satisfied  if 

a*  =WflB(c) 

= 0.  40547  W 

and  from  Eq.  (1 3)  that 


L* -0.20273  W 

In  practice,  we  would  generate  the  random  pulse  height,  , by 
generating  the  random  variable,  > where  has  the  gaussian  proba- 

bility  density 


(I4,  - L)» 

p^  (I^)  = (2no^*H'*  exp  [-i  - 


- 8 - 


and  generate  from  our  random  choice  of  by  the  relationship 

1^,  = i;  exp  (LJ  (20) 

With  this  method  of  choosing  the  random  pulse  height,  ^ , the  mean  pulse 

height  will  be  independent  of  pulse  width,  while  the  second  moment  of  pulse 
height  will  be  so  dependent  on  pulse  width  that,  in  combination  with  the  depen- 
dence of  mean  pulse  count  on  pulse  width  [cf.  Eq.  (4)],  the  k'*'  nature  of 
the  image  power -spectrum  is  insured. 

At  this  point,  having  completed  our  definition  of  a method  of  gener- 
ating the  random  sample  image,  we  need  to  prove  that  this  method  does 
indeed  yield  an  image  with  a k**  power- spectrum.  We  take  this  up  in  the 
next  section. 
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lm>ge  Power  Spectrum 

The  image  generated  by  the  preceding  process  is  a random  selection 
from  an  ensemble  of  possible  images,  the  nature  of  the  ensemble  being  such 
that  the  image  power  spectrum  presumably  behaves  as  . In  this  section 
we  shall  show  that  it  does  indeed  behave  in  this  manner. 


The  randomly  selected  image  can  be  written  as 


Ns 

sft.)  = T 1 9 (piW  . t,_,)  , (21) 

NaO  HaO 

where  we  have  used  tiie  dual  subscript  on  I^  ^ to  indicate  a separate  choice 
of  I for  each  pulse  of  width  a^  , the  choice  being  random  over  the-set 
M s 0 to  . We  let  p = 0,l,2,..,,  1023  denote  the  ppsitions 
in  the  field -of -view.  The  quantity  is  a Poisson  distributed  random 
variable  with  mean  value  1^^^  , calculated  in  accordance  with  Eq. 's  (4) 

and  (7).  The  quantity  is  a log-normally  distributed  random  variable, 

calculated  from  Eq.  (20)  with  L a gaussian  random  variable  selected  in 
accordance  with  Eq.  's  (17),  (18),  and  (19).  The  quantity  is  a gaus- 

sian distributed  random  variable  selected  in  accordance  with  Eq.  (1).  The 
function  P(p;W,  f)  denotes  a pulse  of  width,  W , centered  at  I . Thus, 


1 “ |p-f|  i 

P(p;W,  = < 

0 It  |p./|  > *1,. 


(22) 


To  prove  that  the  form  of  S(p)  has  a power  spectrum,  we 
start  by  considering  the  fourier  transform  of  ^ . Thus  we  write 


- 10  - 


^(H;W,i)=  J dp  exp  (-inp)  ^(p:W.  jl) 


i+a^j/2 


= J dp  exp  (-inp) 


l-^JZ 


= exp  (-inX)  |*  dp'  exp  (-iHp') 


= exp  (-in  X)  a 


sin  (h  a,^/2) 
« H a^/2 


which  is  a random  variable  inasmuch  as  X is  a random  variable. 

Based  on  this,  we  can  write  the  fourier  transform  of  the  randomly 
generated  image  signal,  S(p)  , as 


S(h)  = J dp  exp  (-inp)  S(p) 


" Z Z J*  ®*P 


NaO  MaO 


7 Nk 


" I I 


NaO  MaQ 


^ sin  (h  a,g/2)  . 

" Z *“  iZ  ®*P 
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Here  the  quantities  in  the  curly  brackets  have  been  isolated  because  they 
contain  all  of  the  random  aspects  of  S(h)  . At  this  point,  to  develop  a 
relationship  between  the  random  fourier  transform  and  the  power  spectral 
density,  we  make  use  of  the  fact  that  the  power  spectral  density  function 
associated  with  the  random  image  signal  S , namely  PSD(k)  can  be 
written  as 

PSD(k)  = 2tt  J dn'  <S*{h)  S{h')  (25) 

From  this,  it  follows  that 

f f sin  (h  a /2)  sin  (k'  a /2) 

PSDW=2wY  

NaO  M^e 

Mm  Nm 

I J)  ) . (26) 

Mae  IKaO 

Here  we  have  separated  the  ensemble  average  over  ^ and  over  ^ 
based  on  the  argument  that  the  intensity  and  position  of  each  element  are 
statistically  independent.  The  outermost  ensemble  average  is  over  the 
choice  of  , the  random  number  of  elements.  To  evaluate  the  ensemble 
average  over  in  £q.  (26),  we  first  of  all  note  that  we  are  not  inter- 
ested in  spatial  frequencies  lower  than  (2  a^)-!'  , i.  e. , in  values  of  k any 

smaller  than 

= 2TT/(2iy,)  = tt/^  (27) 

We  shall  evaluate  the  ensemble  average  over  i by  separately  treating 
the  two  cases  in  which  (W,M)  = (W',M'),  and  when  it  is  not.  We  argue  that 
the  case  when  (W,M)  makes  a negligible  contribution,  since  in 

this  case  and  ^ are  statistically  independent  gaussian  random 

variables,  so  that 


- 12  - 


<exp  [i(K 


i 


li 


= exp  [-^(»t3  + K'8)  a,3  + i(K-H')T] 


(28) 


Considering  that  o,  =2  , it  follows  that  for  all  values  of  k of  interest 

to  us,  is  greater  than  2tt*  , and  then,  since  exp  (-2tt*)=  2.  67  XlO"®  , 

the  value  of  the  ensemble  average  is  negligibly  small. 

This  leaves  us  with  only  the  case  where  (W,M)  =(W',M')  as  a possible 
contribution  to  the  power  spectral  density.  In  this  case,  since  and  ^ 

are  identical,  we  can  write 

{(exp  [i  (h  - K u»)]  ) ]k. K#, 

= exp  [ (H-K')a  a*  + i (h  -h')  T]  (29) 

Because  of  the  relatively  large  value  of  a,  , we  can  treat  the  gaussian  in 
Eq.  (29)  as  a delta  function,  i.  e. , as  (2n  6(k-h')  in  evaluating  the 

h'- integration  in  Eq.  (26).  Thus  we  get 

7 sin  (h  \/2)  ^ 

PSD<k)  = (2n)«  c,  I K)*[-77;72 • P”* 

M>0  MaO 

where  here  the  ensemble  average  is  only  over  the  choice  of  for  each 
value  of  W , and  the  choice  of  ^ for  each  value  of  (W,M)  . These 
two  random  features  are  statistically  independent,  and  factorable  in  terms 
of  their  effect  on  the  power  spectral  density  as  defined  in  Eq.  (30).  Thus, 
since 

<V«*>=<V>  (31) 
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1 


i.  e.  , the  value  of  is  independent  of  M , (though  as  we  have  set 

up  the  statistics  of  I , it  retains  a dependence  on  W ),  we  can  rewrite 
Eq.  (30)  as 


PSDt.)  = (2r)«  ^ (..)•  <!,•>  [-,  J 


(32) 


MaO 


Making  use  of  £q.  's  (3),  (4),  (7),  (8),  and  (9).  and  the  statement  immed- 
iately following  £q.  (12),  we  can  rewrite  £q.  (32)  as 


PSD(k)  = (2TT)a'a(,,  \ a^*  2-^  b**  c-“ 

N*0 

, -<«♦!) 
sin  (h  a^  2 ) . 

H ^ 2-(-*n  J 


_ ’ sin  (h  a.  2-( "♦!>)  • 

(2")«  c,  [ 2-  [— TT^P^] 


HaO 


(33) 


It  is  easy  to  see  that  the  power  spectral  density  is  the  sum*of  a set  of 
seven  power  spectra  each  of  which  is  constant  at  low  frequencies,  but  then 
breaks  at  some  "knee-frequency",  and  after  that  falls  off  at  12  dB  per 
octave.  Each  of  the  seven  successive  spectra  has  a constant  region  value 
that  is  a factor  of  two  less  than  the  preceding,  but  a "knee -frequency" 
that  is  a factor  of  two  higher.  Since  12  dB  per  octave  corresponds  to  a 
factor  of  four  fall-off  for  each  factor  of  two  increase  in  frequency,  it  is 
apparent  that  at  each  knee  frequency  there  is  some  one  of  the  sequence 
of  seven  subsidiary  power  spectra  that  is  dominant,  with  the  one  just  pre- 
ceding and  the  one  just  following  reduced  by  a factor  of  two.  Thus  the 
stun  of  all  the  subsidiary  power  spectra  should  vary  with  frequency  as 
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the  response  at  the  knee -frequency  varies  with  knee -frequency.  This  is 
• K**  dependence.  This  situation  is  indicated  in  Fig.  1,  where  we  show 
the  curves  for  W = 0 , 1 , 2 , . . . , 6 (W  = 7 is  off  the  graph),  and 
the  sum  of  these  curves  making  up  the  power  spectral  density.  The  sum 
quite  clearly  follows  the  knees  of  the  W-curves,  and  therefore  has  the 


hH.  (Jepe  ndenc  e . 


5. 


CommentB 


With  the  technique  defined  in  the  preceding  section,  we  have  estab- 
lished a method  of  generating  a one -dimensional  image  with  a power 

spectrum.  Our  next  objective  will  be  to  use  this  procedure  to  develop  sample 
images  which  we  will  then  filter  to  simulate  the  effect  of  the  CIS  transfer 
function,  and  after  that  corrupt  with  Poisson  distributed  noise,  i. e.  , re- 
placing the  filtered  intensity  at  each  signal  point  with  a Poisson  distributed 
random  variable  with  that  mean  value.  We  will  then  subject  such  a one- 
dimensional image  to  the  nonstationary  post-detection  compensation  pro- 
cedure which  we  described  earlier. 
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exact  ^sin3(x)/x]^  dependence  for  , but  in  the  oscillatory  region  we  have  suppressed 

this  oscillation  and  replaced  the  curve  with  the  dashed  line,  which  is  its  envelope.  The 
envelope  curves  are  shown  for  , 1 , ....  6 . The  dot-dash  curve  is  the  sum  of 

these  envelope  curves  and  represents  the  power  spectral  density.  Its  dependence  is 
apparent  for  spatial  frequencies  above  h 1/a  . 


mm'. 
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